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ABSTRACT 



G> ' Assuming a diamagnetic interaction between a stellar-spot originated localized mag- 

netic field and gas blobs in the accretion disk around a T- Tauri star, we show the 
possibility of ejection of such blobs out of the disk plane. Choosing the interaction ra- 
^ I dius and the magnetic field parameters in a suitable way gives rise to closed orbits for 

the ejected blobs. A stream of matter composed of such blobs, ejected on one side of 
the disk and impacting on the other, can form a hot spot at a fixed position on the disk 
^ ! (in the frame rotating with the star). Such a hot spot, spread somewhat by disk shear 

[^^ I before cooling, may be responsible in some cases for the lightcurve variations observed in 

0> ■ various T- Tauri stars over the years. An eclipse-based mechanism due to stellar obscu- 

^ . ration of the spot is proposed. Assuming high disk inclination angles it is able to explain 

Y^'. many of the puzzling properties of these variations. By varying the field parameters and 

O I blob initial conditions we obtain variations in the apparent angular velocity of the hot 

-*-» ' spot, producing a constantly changing period or intermittent periodicity disappearance 

ci • in the models. 

>'■ 

^ ■ Subject headings: stars:formation — stars: TTauri — accretion, accretion disks — mag- 

^ ! netic fields 
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I. Introduction 

Most T-Tauri stars (TTS) exhibit irregularities 
in their luminosity, observable photometrically and 
spectroscopically. This variability is observed on a 
variety of timescales and in a wide frequency band. 
Although the short-term variability can be explained 
as variable chromospheric eruption/emission episodes 
(Wordcn et al. 1981), the long-term variability is not 
yet fully explained. An extreme example is BP Tau, 
whose lightcurve, monitored over 50 days (Simon et 
al., 1990), shows periodic brightness fluctuations with 
a period of 7.6d, then an almost constant brightness 
level for 2 weeks, and finally the periodic behavior 
resumes - but with a period of 6. Id. Variability of 
this type is well documented for a variety of stars 
(e.g. Herbig 1962, Bouvier et al. 1995, Rydgren & 
Vrba 1983). The puzzling property of a number of the 
classical TTS (CTTS) is that their variability period 
changes on a time scale of weeks to months (BP-Tau, 
mentioned above, is the only case where the period 
change was actually continuosly monitored). 

One type of explanations of this effect involves the 
stellar magnetic field, following the idea of Ghosh & 
Lamb (1978) (see also Konigl, 1991). An oscillation 
of sorts is assumed to occur in the Alfven radius and 
thus the magnetic boundary layer there changes its 
size and position, causing the variability (Smith et 
al. 1995, Armitage 1995). The weakness of these 
models lies in the need for an organized strong dipole 
magnetic field in TTS, which is rather doubtful obser- 
vationally (Menard 1996). Another promising model 
is that of cool and hot spots on the stellar surface, 
with cool spots analogous to sunspots and hot spots 
being the by-product of the accretion process (Vrba 
et al. 1986, Bouvier & Bertout 1989, Bouvier et al. 
1993, 1995). 

A recent series of observations, the COYOTESQ 

II, were obtained over a 2-month monitoring of the 
light variations of TTS of the Taurus-Auriga dark 
cloud. The results provide further evidence for tem- 
poral variations of the photometric periods of CTTS. 
For example, IQ Tau and GM Aur exhibit relatively 
flat lightcurves with narrow minima, while the other 
lightcurves are more sinusoidal. DR Tau is a special 
case, as both the length and depth of the minima vary, 
and the lightcurve is neither sinusoidal nor flattened, 
but exhibits mixed traits. In the CW Tau case, the 
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star becomes bluer in U-V and B-V when fainter than 
V=14.0'". It also exhibits large photometric variabil- 
ity, with amplitude variations ranging from almost 3™ 
in the U-band to 1.5™ in the I-band. 

Several photometric periods were reported for DF 
Tau at different epochs: 8.5d (Bouvier & Bertout 
1989), 7.9d (Richter et al. 1992), and 9.8d in the 
COYOTES II observations. DR Tau and DG Tau 
may be also candidates for this category. In the COY- 
OTES I campaign two possible periods were found 
for DR Tau (2.8d and 7.3d), while the COYOTES II 
data points to 7.23d; the DG Tau data from COY- 
OTES I/II reveals a possible period of 6.3d, while 
a period of 4.3d is claimed by Guenther(1994, priv. 
comm with J. Bouvier). From the data for TAP 57NW 
it also seems that the period and the overall shape of 
the lightcurve vary with time. In addition to these 
period variations, the same stars sometimes do not 
exhibit any periods in their lightcurve (e.g. DF Tau 
- Rydgren et al. 1984, Bouvier & Bertout 1989, BP 
Tau- Simon et al. 1990), indicating that the period 
not only changes but may disappear temporarily. 

Only stars whose lightcurves are consistent with 
the presence of hot spots were found so far to ex- 
hibit period changes and disappearances. This sug- 
gests that the cause of these phenomena is to be found 
in the accretion process itself. This is also suggested 
by the spectroscopic signature of mass inflow at free 
fall velocities (Hartmann et al. 1994, Edwards et al. 
1987), which provide support for the magnetic accre- 
tion mechanism. 

The period variations cannot be attributed to true 
changes in the stellar rotation, as the response time of 
the stellar angular velocity to non-steady accretion is 
of the order of lO^y (Pudritz & Patel 1994) while for 
example the timescale over which the period changes 
are observed to occur is two weeks for BP Tau (Simon 
et al. 1990). Thus one is bound to conclude that the 
cause of the period variations lies in changes of the 
angular velocity of the hot spots themselves. 

2. The model 

King & Regev (1994) (hereafter KR) proposed that 
the low rotation rates and outflows in TTS could be 
explained given the presence of a magnetic loop struc- 
ture on the central star in a T Tauri-disk system. The 
loop is assumed to interact with the disk material via 
a magnetic surface drag force exerted on diamagnetic 
gas blobs, following the scheme of Drell, FoUey and 



Ruderman (1965). These authors showed that the 
energy-loss timcscale for a diamagnetic object of mass 
m crossing the lines of a magnetic field B is 
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with I a typical length scale for the object and ca the 
Alfvcn velocity in the plasma outside the object. The 
minimum conductivity condition for this equation to 
hold is 

^ > 7T-T- • (2) 



2TrlcA 

This reduces to a condition on the blob Icngthscale 
which is satisfied for all realistic parameters appro- 
praite for TTS conditions (see KR). 

Defining now a drag coefficient k ~ t'j^ran' ^^*^ drag 
force per unit blob mass is written as 



fd™s--fc[v'-(v'.b)b], 



(3) 



where v' is the blob velocity relative to the field lines 
and b is a unit vector in the direction of the field, 
fdrag IS thus perpendicular to the field lines and its 
manitude is proportional to the perpendicular veloc- 
ity component. The drag coefficient, k, depends on 
the magnetic field strength and on the local interblob 
plasma density (through ca)- Therefore it is a func- 
tion of position. In addition, it depends also on the 
individual blob parameters. 

We shall use a magnetic field of a flux tube re- 
sulting from a stellar magnetic loop crossing the disk, 
with the following functional dependence on position: 



B{x,y,z) = BqC 
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where x, y, z are Cartesian coordinates (the origin is 
at the star's center and z = is the disk plane), Bq 
is a constant and Wp , Wz are the Gaussian spreads 
characterizing the loop. 

The functional form of the field should approxi- 
mate the spatial dependence of a magnetic field lo- 
calized in a flux tube with a horizontal spread Wp 
around the point a;o, j/o- A Gaussian form is chosen as 
a convenient possibility. The vertical cutoff (Gaussian 
decline in the field strength, with spread Wz) reflects 
our wish to avoid repeated blob-flux tube interaction 
in our simple model. In any case, a (Gaussian) ver- 
tical cutoff in k{r) will follow from its dependence on 
the disk density (see below) , which is usually assumed 
to have a Gaussian decay with height. 



The full expression for the drag coefficient, k, is 
obtained by substituting ca = B/^/4jTp and m = Ipb 
in the original definition, pb and p are the densities 
of the blob and the interblob plasma in the disk, re- 
spectively. Thus 



k = 2y^Bp^ —- . 
Pbl 



(5) 



The equations of motion of a blob are written in a 
frame rotating with the star, with angular velocity uJ 
as 

GMr „^ . ^ ,^ ^ , /. .. ,-x,'\ 
r~ g 2a; X r — (jj X (w X r) — fc I r — (r • b)bj , 

(6) 
with r = (x, y, z) and b = (bx, by, bz) defined as 

6a; =: sin cos , fej, = sin sin ,bz ^ cos 9 . (7) 

9 and are angular parameters corresponding to the 
flux tube direction relative to the disk. M is the cental 
star's mass and G the gravitational constant. 

We first assume that lu = ujz, i.e. the rotation axis 
of the star is the z-axis of the coordinate frame (z is 
the corresponding unit vector). This is very reason- 
able and means that the star and disk rotate around 
the same axis. Next we nondimensionalize the equa- 
tion of motion, scaling the physical variables by the 
following chractcristic units: distances are scaled by 
the value of the corotation radius, defined here to be 
the radius at which the local Keplerian angular ve- 
locity equals Lu-the angular velocity of the star and 
thus the field, rco = {GM/uj'^y^^; the time by the ro- 
tational time T — l/uj and the rotational velocity by 
LU. The following nondimcnsional equation of motion 
then results: 



r = T^2z X r — z X (z x r) — aK{r) ( r — (r • b)b) , 

(8) 
where the nondimcnsional drag coefficient is 
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Here we have assumed that the interblob disk density 
depends on z only and po is the midplane value. As 
mentioned before this causes the appropriate decline 
of K with height. The constant parameter a is 
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where P* = 27r/a; is the rotation period. With the 
various physical variables expressed in their typical 
values the parameter a can be written as 



a = 0.4875 



Po 



lOOG V 10^10 [5/cto3] 

I Pb 



1/2 



P* [days] 



109[cm] lO-^g/cm^] 



(11) 



The interaction is thus parametrized by the inter- 
action strength a, the flux tube width Wp and the 
tube direction 9, (f>. Wz should be approximately the 
disk half-thickness. The motion of a blob after suffer- 
ing an interaction with a magnetic loop will obviously 
depend also on the initial blob velocity relative to the 
field line. As a rule the blobs will be taken to orbit 
the central star with a local Keplerian velocity plus a 
small inward drift, hence the importance of the inter- 
action radius tq = -^/xq + j/q. 

Eq. ^ describes the blob motion in a rotating 
frame, with the Coriolis and centrifugal terms in- 
cluded as appropriate. It can be solved numerically 
for a variety of initial conditions and loop and blob 
parameters. We shall refer to the solutions of this 
equation as ballistic orbits if the trajectory leaves the 
disk and crosses its plane again. 

3. Ballistic orbit calculations 
3.1. General considerations 

We start the integration of the blob motion from a 
point close to the interaction point (xq, j/o) in the disk 
plane, but far enough from it so that the magnetic 
drag term is negligible. A distance of several Wp is 
appropriate and thus with Wp = erg and e <C 1, we 
obtain a sharply (with respect to tq) increasing field. 

For ro < Tco the blobs are slowed down by the mag- 
netic drag (see KR). If no vertical velocity component 
is acquired this will merely increase the accretion rate. 
However, significant vertical velocity component may 
often be produced by the interaction; in such cases 
the blob leaves the disk - only to fall back on it (close 
to the star in most cases), after a very inclined orbit. 

In our calculations the blob is started with a Ke- 
plerian velocity (in the inertial frame) plus a small 
inward drift. For the case ro < Tco the blob is started 
behind of flux tube, so as to let the blob overtake the 
tube. In this work we shall concentrate only on this 
case, since rg > Tco implies almost always ejection of 



the blob out of the disk plane without subsequent im- 
pact, see KR and Pearson & King (1996). As pointed 
out in the Introduction, we focus here on the possi- 
bility of blob ejection out of the disk plane, but with 
less than escape energy and with the blob supposed 
to land quite close to the star (see below). Indeed 
we were able to obtain such a situation, for a large 
fraction of the parameter space, in the case rg < Vco- 

For a given choice of the parameters, each ejected 
blob will follow exactly the same path and thus its 
center of mass will impact the disk at the same posi- 
tion (in the rotating frame) as all the preceding ones. 
If we envisage a succession of blobs being ejected one 
after the other and forming during their flight a con- 
tinuous stream of matter, this stream will impact the 
disk at a fixed (in the rotating frame) position. The 
impact point with the aforementioned hot spot cre- 
ated by the stream will thus move on the disk, around 
the star, with an angular velocity to. If the spot is 
close enough to the star for a given inclination angle, 
it will be, at least partly, eclipsed by the star. The 
observed luminosity variability due to this spot will 
then have the fixed frequency lj. The stream may also 
impact the star and the observational result, regard- 
ing time variability, should be the same. Now, if we 
allow the field parameters to vary in time, the spot is 
bound to move, in the rotating frame, on a timescale 
of the order of the magnetic field variation timescale. 
We shall assume the this timescale approximately of 
the order Ib ^ 5P*, where P* = 2t: /lu is a rotation 
period of the field. Thus we postulate that the field 
configuration changes (for example due to the inter- 
action with the disk plasma and the blobs, or due to 
internal changes) , and this happens on a time scale of 
~ 5 rotations. A spot moving in the rotating frame 
can explain all the unusual features in the observed 
light curves. The period can change on this time scale. 
The variability may disappear (if the spot moves far 
enough out as not be eclipsed any more) and reappear 
again, if the spot moves back in. A systematic study 
of all the possibilities calls for a many body particle 
numerical simulation for the blobs in the disk. In ad- 
dition we need to simulate numerically the stream of 
the ejected matter and its impact on the disk. Such 
a project is now in progress, while here we would like 
only to demonstrate the idea by describing various 
cases of individual blobs. 

If we fix all the parameters but start the calculation 
with different initial conditions (due to blobs coming 
from slightly different places), we expect a spread in 



the impact point in the rotating frame, giving rise to 
an extended hot-spot. The typical spread of impact 
points we have found is of the order of the stellar 
radius, i?* for the cases described below, with the 
initial blob conditions around the magnetic flux tube 
set to give a ballistic orbit. 

3.2. Results for varying field parameters 

The canonical values of the parameters for our cal- 
culations are: = 7r/4, cj) = 37r/4, r^ — 0.85, a — 100 
in non-dimensional units, with time scaled by P^ and 
lengths by the value of the co-rotation radius Tco- We 
chose always a single blob initial condition, remem- 
bering that a spread in the spot size is expected due 
to blobs with close initial conditions. We then vary 
(linearly in time) each of these parameters in turn, 
keeping the others fixed. The timescale of change is 
of the order of 5P^. The Gaussian spreads are always 



kept fixed at Wp 



0.05. 



The variation of a has a relatively small effect on 
the impact point position in the rotating frame. Thus 
the spread in blob sizes and masses as well as the 
actual absolute value of the field matter little in this 
model (A variation of a between 1 and 10^ resulted 
in a change of only ~ 2% in spot angular velocity, see 
figure |(a). 

Significant change occured when we have varied 
the value of rp, the radius of the point where the loop 
crosses the disk (in units of Tco)- A linear change 
ro = 0.6-|-0.08t (variation in the range 0.6 — 0.92 dur- 
ing 4 rotations) resulted in changes of the spot center 
position in polar coordinates {Rimpi^imp), shown as 
a function of time in figure Mb). Unlike in the vari- 
able a case, (figure ^a)), here the apparent angu- 
lar velocity of the spot changed drastically; combined 
with the spot grazing the star itself part of the time 
{Ri, ~ 0.23rco) this produced a lightcurve very much 
resembling that of DR Tau (see fig. 3 of Bouvier et 
al. 1995). The model lightcurve exhibits variation in 
shape, period and depth of the minima (figures 0,0) 

By varying cf) and 9 we again obtain significant mo- 
tion of the spot in the rotating frame. Figure n depicts 
these cases as well: spot's coordinates evolution with 
9 varying as 6* = 0.5 -I- 0.2t- curve (c), and the coor- 
dinate evolution with cj) = 1.9 + 0.3t- curve (d). We 
can see that a spot moves away from the star, and so 
the possibility of period disappearance and its conse- 
quential return is present, although with a different 
apparent angular frequency. It should be also noted 



that although fi (spot's angular velocity as seen by an 
observer in an inertial frame) varies some 10 — 15%, 
the angular velocity deduced from the spectrogram is 
just the average. In both cases the apparent period 
differed by about 10% from P-,.(fig. |) 

3.3. Model predictions 

In order to be able to predict whether the vari- 
ability can be observed, one needs to know the stellar 
radius, i?*, the angle of inclination of the disk and the 
spot size. The apparent angular velocity of the spot 
is clearly ft = uj + Qimp- This can be calculated with 
no further assumptions. We assume that the spot size 
is very small {^ 0.2_R^,) and the inclination angle is 
such that eclipses are always occurring for the cases 
studied. The star is assumed to have i?* = 2i?0, 
fli, = uj = O.irJ/^-^ (a tenth of the breakup value). 
The results should thus be considered as illustrative 
only. In figures g,|| we see the angular velocity of the 
spot in the inertial frame, i.e. as should be observed, 
the and the power spectrum of the predicted inten- 
sity variation- for the three cases of varying different 
parameters (corresponding to the results depicted in 
figure 0). 

4. Discussion 

The above model is an extremely simplified view 
of the whole problem. In reality, a process for the 
creation of the diamagnetic blobs would have to be 
specified; the interaction between blobs should be in- 
cluded; and the disk thickness has to be taken into ac- 
count, as it is unreasonable to expect that a blob flung 
out from the downward side upwards will traverse the 
disk without colliding with the disk and other blobs 
inside it. Here we do not specify the MHD insta- 
bility responsible for the creation of blobs, but only 
note that in magnetic Cataclysmic Variables there is 
probably observational evidence for a blobby accreted 
plasma (see Wynn and King 1995). In addition, we 
need the blobs to exist only for a time sufficient for 
the magnetic interaction in the disk. 

If the conditions confining the blobs are specific to 
the disk, they will expand on the thermal time-scale 
after ejection. Expanding, they will merge together 
into a stream-like structure; thus the above picture of 
single blobs falling into the disk is useful only in the 
context of checking the whole idea; the results should 
be taken as no more than qualitative and representa- 
tive of the general scales. 



Still, the model looks promising, as it might ex- 
plain the variability of the observed period in certain 
TTS; for example, the disappearance of the period al- 
together can be explained by the hot spot moving far- 
ther away from the star, thus evading being eclipsed. 
The model works only for highly inclined systems. It 
will be thus interesting to look for an observational 
correlation between period variability and inclination 
in CTTS. We do not rule out stellar spots as causing 
variability in the case when disks are absent (then 
this is the only effect) or when they are present (then 
the star and disk spots are both present, but only the 
disk spots are responsible for period changes). 

In the next step of this research we intend to per- 
form a three dimensional SPH simulation of the accre- 
tion disk and the stream resulting from the magnetic 
interaction taking into account in the mass infall rates 
as well (which shall affect the general spot luminos- 

ity). 

Concluding, we make the following remarks: 

1. The optical thickness of the stream and the pen- 
etration depth in the disk are functions of the 
stream density. We expect that if the stream 
density gets smaller, the observed temperature 
will get closer to the shock temperature; mean- 
while the total luminosity will decline. This 
could explain the observed blue shift during the 
low luminosity phase in some TTS. 

2. In some instances the observed period disap- 
pears altogether, and reemerges afterwards. In 
our model, for low disk inclination angles i, the 
spot will not be eclipsed for Rimp > R*/ sm{i). 
As we see in the examples, Rimp varies in a wide 
range, so this is quite plausible. Another possi- 
ble explanation is that the stream disappears al- 
together (for example if the loop lies in the disk 
plane), and reappears on the other side of the 
disk, its luminosity greatly diminished by the 
disk thickness. In addition, although improb- 
able, the spot angular velocity could be very 
close to if a suitable parameter variation oc- 
curs. The spot could thus be visible for a long 
duration. 

3. We do not take here into account the scenario 
where the spot falls on the star itself. This oc- 
currence is quite possible, especially when tak- 
ing into account the spot size. At least a part 
of the spot may be on the star, while the other 



part is still on the disk. For instance, a disk spot 
changing into a star spot will have the follow- 
ing influences on the lightcurve: as a growing 
proportion of the stream falls on the star, the 
apparent area of the spot will grow; thus the lu- 
minosity will rise. Being on the stellar surface, 
the spot apparent size will vary as sm{il,spott), 
and the luminosity will vary sinusoidally. We 
may expect that the lightcurve will look some- 
what similar to the one we obtained in figure 

i 

4. If the spot grows beyond 2i?*, we shall en- 
counter shallower minima as only a part of the 
spot will be eclipsed. Taking into account the 
fact that a star is spherical, even for a constant 
spot size the eclipsed part shall decline as the 
spot falls further from the star, again resulting 
in shallower (and shorter) eclipses. 
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Fig. 1. — Spot coordinates (in the rotating frame) as 
a function of time (a) for varying a [a in the range 
1-10^), (b) ro(range 0.6-0.92 rco),(c) </> (range 1.9-3.1 
rad),(d) Q (in the range 0.5-1.3 rad). t - time in P^. 
units 

Fig. 2. — Apparent spot angular velocity, with the 
same parameters as figure 1. We can see that the 
maximal variation occurs in the tq case (b), while 
varying ot hardly infuences the results (a). Variation 
of (c) and d (d) produce a comparable variability of 
some 10-15%. (t - time in P^ units) 

Fig. 3. — Power spectrum of the calculated intensity 
variations. (a)In the varying a case the frequency 
obtained equals 0.9 27r/ri^. (b) Varying Rq case. The 
deduced frequency is half the actual one, as one of the 
minima would be inobservable. (c) Varying (\) case. 
The power spectrum suggests a frequency that is 10 % 
less than the actual, (d) Varying Q case. The deduced 
frequency is about 15% more than the actual, (t - time 
in Pt units) 

Fig. 4. — Modelled spot luminosity (variable tq case). 
Both the period and the eclipse length vary, as does 
the overall shape of the light curve. The second min- 
imum (arrow) is too shallow to influence the power 
spectrum. 



This 2-colurnn preprint was prepared with the AAS lATJiiX 
macros v4.0. 
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